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Abstract 
Summary: 

We describe libRoadRunner, a cross-platform, open-source, high performance 
CH — h library for running and analyzing SBML-compliant models. libRoadRun- 
ner was created primarily to achieve high performance, ease of use, portability 
and an extensible architecture. libRoadRunner includes a comprehensive API, 
Plugin support, Python scripting and additional functionality such as stoichio- 
metric and metabolic control analysis. 

Accessibility and Implementation: 

To maximize collaboration, we made libRoadRunner open source and released 
it under the Apache License, Version 2.0. To facilitate reuse, we have devel- 
oped comprehensive Python bindings using SWIG (swig.org) and a C API. Li- 
bRoadRunner uses a number of statically linked third party libraries including: 
LLVM [4], libSBML [1], CVODE, NLEQ2, LAPACK and Poco. LibRoadRunner 
is supported on Windows, Mac OS X and Linux. 

Contact: 

hsauroOu . Washington . edu 

*to whom correspondence should be addressed 
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Supplementary information: 

Online documentation, build instructions and git source repository are available 
at http://www.libroadrunner.org 

1 Introduction 

Over last few decades researchers have developed many tools to simulate bio- 
chemical networks [7] . Most of such simulators are however embedded in mono- 
lithic environments, making it very difficult to reuse the code in other appli- 
cations. To resolve this issue various groups have built simulators as reusable 
libraries. Most notable of these are COPASI [2], libSBMLSim [9], RoadRunner 
C# [7], SOSLib [5] and SBSCL [3]. Despite all these efforts, finding a simulator 
that is easy-to-use, has good SBML-compliance, is capable of handling complex 
models and fast enough to run thousands of models is challenging; for example, 
in multi-cell simulators such as CompuCell3D. We designed libRoadRunner to 
address these shortcomings and by using Just-In-Time (JIT) compilation via 
Low Level Virtual Machine (LLVM), we have significantly improved computa- 
tional performance as compared to existing solutions. 

This version of libRoadRunner has a new internal design compared with the 
previously developed C# roadRunner. The APIs are designed with the modeler 
in mind, offer a wide range of functionality and provide high performance access 
to model and simulation data, making libRoadRunner ideal for realtime inter- 
active environments. To facilitate easy extensibility libRoadRunner provides a 
plugin API making it possible to easily write extensions. Figure 1 shows the 
architectural overview of libRoadRunner. 

The most novel and important component of the libRoadRunner is its com- 
putational core based on the LLVM. LLVM is designed for real-time optimiza- 
tion and dynamic compilation of computer languages. This has allowed us to 
achieve extremely low runtimes. This speed is particularly important for large 
simulation applications such as multi-cell models that require many thousands 
of biochemical networks to be simulated simultaneously. libRoadRunner 1.0 
supports all models in the SBML test suite except for those that include delay 
equations and algebraic rules. 

2 Methods 

2.1 LLVM Backend 

We have developed a non-tracing, multi-pass JIT compiler for the SBML lan- 
guage using the LLVM framework which allows us to directly JIT-compile SBML 
in memory into native machine code at runtime. Our compiler treats SBML as 
a statically typed, dynamically scoped declarative language allowing us to not 
only generate specific code for each SBML model, but to also generate a specific 
data layout that is optimized for each model. As a result the model and the 
integrator share the same data layout eliminating the need to copy data between 



2 



Downloaded from http://biorxiv.org/on September 18, 2014 



the integrator and the model evaluation code. The flexible architecture of the 
compiler also greatly facilitates additions of SBML extensions 

2.2 Extensibility 

A key design requirement of libRoadRunner was to provide a component-based 
architecture. Therefore, internal components such as the integrator, steady 
state solver, SBML compiler, etc., only communicate with each other through 
pure virtual interfaces. This design allows us to readily add new internal com- 
ponents such as a planned GPU SBML compiler, new integrators, stochastic 
integrators and steady state solvers. While the easiest way to extend function- 
ality is to write Python scripts using the Python API, for performance sensitive 
applications developers can access the same functionality through the C++ or 
C API's. We also offer a C Plugin API that allows developers to write dynam- 
ically loadable plugins in any C complaint language such as Fortran or Object 
Pascal. The Plugin API is written such that it allows generic access to all 
plugins. For example, scripting languages such as Python can access a plugin 
without having to implement a new Python API for every plugin. The current li- 
bRoadRunner distribution comes with two example plugins that illustrate how 
libRoadRunner functionality can be extended. One example implements the 
Levenberg-Marquardt algorithm [6] for fitting models to experimental data. 

2.3 Python API 

LibRoadrunner provides a comprehensive and pythonic Python module which 
is designed to have the feel of the Python SciPy package (scipy.org) and exposes 
the full functionality of the public C++ API. Consequently, the easiest way 
to use libRoadRunner is via Python scripting. This API allows access to any 
SBML element along with calculated values such as the Jacobian, eigenvalues 
or the various control coefficients using the Python dictionary protocol. All 
libRoadRunner data arrays are accessed directly without copying as Numpy 
arrays. This provides a commonly used access protocol. Using the Python API, 
modelers can create one or many RoadRunner simulator objects which offers 
a very natural way of using libRoadRunner as the Reaction-Kinetics engine in 
multi-cell simulators such as CompuCell3D (www.compucell3d.org). 

2.4 Performance 

To demonstrate the capabilities of libRoadRuner we compared it to three simu- 
lator libraries: libSBMLSim, COPASI and SBSCL. We tested the ability to pass 
the SBML test suite (1008 models), the ability to simulate very large models 
and the ability to simulate systems with large numbers of events. Table 1 lists 
the times recorded from four different libraries based on three types of models. 
The first model tests the ability to scale to large numbers of state variables 
defined by relatively simple ODEs. The second model tests the ability to deal 
with complex SBML functions. The timings are the total wall time to complete 
each process. Tests were run on a 2.6 GHz Mac Pro, OS X 10.6, and used 
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Figure 1: Architectural Overview of libRoadRunner. 

COPASI v. 4.9.43, LibSMLSim v. 1.1 and SBSCL v. 1.3. Note, we experienced 
instability of LibSBMLSim runs with >100 copies of the Brusselator system or 
complex SBML. 



2.5 Using libRoadRunner in multi-cell simulations 

We have successfully deployed libRoadRunner as the Reaction-Kinetic solver en- 
gine in CompuCell3D - a simulation environment used to build multi-cell based 
models of tissues. We have used the libRoadRunner Python API to associate 
each cell with a collection of libRoadRunner-based simulators. Since Compu- 
Cell3D uses Python scripting to describe its models, integrating libRoadRunner 
was straightforward. 



Table 1: All times are in seconds. The first set of tests consisted of N copies 
of the Brusselator system in a single model. The second set was a sin func- 
tion implemented as a 63 element piecewise SBML functions combined with N 
parameters defined by rate rules integrating the sin function. Test models are 
available at libroadrunner.org. 
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We ran the Delta-Notch patterning simulation [8] on an array of 36 adjacent 
cells with different reaction-kinetics solvers and determined that the libRoad- 
Runner solver led to a two fold improvement in runtime. In multi-cell simula- 
tions time spent on solving reaction-kinetics models makes up only a portion of 
the total runtime. Thus the two fold speedup in the overall simulation time is 
significant. 

3 Discussion 

Decades of development has produced powerful, intuitive and feature-rich soft- 
ware to simulate biochemical networks. The existing portfolio of SBML-compliant 
tools satisfies the needs of most modelers. Yet, limitations of existing packages 
present a significant research barrier for applications that require high perfor- 
mance, excellent SBML compatibility, broad functionality and extensibility. Li- 
bRoadRunner addresses these shortcomings. It is designed for researchers who 
need a simulation engine for interactive desktop simulation, multicellular sim- 
ulations, which often require 10,000 or more simultaneous simulations, and for 
web based interactive simulation with multi-session requirements. Even though 
the current version of the RoadRunner appears to be complete, we are plan- 
ning a number of additional features including even better SBML compliance, 
stochastic GPU integrators, parallel distributed model JIT compilation and sim- 
ulation, spatial models as well as new plugins such as alternative optimizers, and 
bifurcation analysis. By using an open-source, open-access development process 
we hope to engage external developers which should result in higher quality and 
more frequent future releases of libRoadRunncr. 
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